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Abstract 

We exhibit a class of integer spin systems whose free energy can be written in term of 
an absolutely convergent series at any temperature. This class includes spin systems on 7L d 
interacting through infinite range pair potential polynomially decaying at large distances r at 
, a rate l/r d+e with s > 0. It also contains the Blume-Emery-Grimths model in the disordered 

phase at large values of the crystal field. 

1 The class of spin systems: definitions and results 



Let V be a countable set. We define an integer spin system on V by supposing that in each site 
x € V there is a random variable a x (the spin at x) taking values in the set {0, ±1, ±2, . . . , ±iV} 
where N is an integer. For U C V, a spin configuration ajj in U is a function a\j : U — > 
{0, dbl, ±2, . . . , ±N} : x t— > a x . We denote by £[/ the set of all spin configurations in U. Given 
a spin configuration a\ in the finite "volume" A C V, the Hamiltonian of the system at volume 
A (with free boundary conditions) is defined as 

o 

ON 

■ {x,y}cA xeA 



H a (oj)= V(x,y,a x ,a y ) + D^2a 2 x (1.1) 



Typically, the set V is the unit cubic lattice Z d , but theorem 1 below holds for any countable 
set V regardless of its structure. The assumptions on the pair potential V(x,y, a x , a y ) in the 



' Hamiltonian (1.1) are the following. 

A. For any pairs x, y such that a x a y = we have 

V(x,y,a x ,a y ) = (1.2) 

B. There exists a positive function J(x,y) and a positive number J such that 
\V(x,y,a x ,a y )\ < J(x,y), V ({x, y}, a x , a y ) (1.3) 

and 

sup V J(x,y) = 2J (1.4) 

Note that assumption (1.4) includes infinite range interactions polynomially decaying in a 
summable way. Note also that, when D > J, the lowest-energy state (i.e. the spin configu- 
ration a A € Sa for which H\(a\) attains its minimum) is the spin configuration a A = 0, i.e. all 
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sites have spin equal to zero. The region D > J is usually called the disordered phase and the 
parameter D is called the crystal field. 

The partition function of the system in the volume A, at inverse temperature (3 is given by 

Z A {(3)= e"^ (CTA) (1.5) 

and the free energy at finite volume A is given by 

fA((3) = -^MZ A ((3) (1.6) 

where here and elsewhere in the paper |A| is the number of elements of the finite set A. 
Our results on the system above can be summarized by the following theorem. 

Theorem 1. Consider the spin system with Hamiltonian (1.1). Under the assumptions (1.2)- 
(l-4)> there exist positive numbers D c , Pi and fa with D c > J and fa < fa such that the free 
energy /a(/3) defined in (1.6) can be written in terms of an absolutely convergent series uniformly 
bounded in A either if (3 € [0,/?i] or (3 S [fa,oo) and D > J. Moreover if D > D c then f\{fi) 
converges absolutely for all (3 G [0, oo). 

We remind the reader that the standard polymer expansion has been shown to be absolutely 
convergent at sufficiently high temperatures or low density for wide class of spin systems, either 
using Kirkwood-Salzburg equations [8] or by bounding directly the coefficients of the expansion 
|17j . Moreover it has been shown in [9] that analyticity of a similar class of spin systems at high 
temperatures or low densities can be inferred from the Dobrushin uniqueness theorem. 

On the other hand it is well known that for a wide class of spin systems a unique zero- 
temperature ground state implies a unique Gibbs state at sufficiently low temperatures and 
hence absence of phase transitions (see e.g. |12| and references therein). Analyticity at low- 
temperature via Pirogov-Sinai theory has been established in [19] for spin systems interacting 
via finite-range potentials and has been extended in [13] to infinite-range potentials. 

Here we present a class of spin systems, depending on a positive parameter D, whose free 
energy is analytic simultaneously at high and low temperatures if D > J or at any temperature, 
if D > D c . It is worth remarking that our class of spin systems also includes the Blume-Emery- 
Griffiths (BEG) model @] on 7L d . To recover the BEG model just choose V = Z d , N = 1 
and 

V(x,y,a x ,(T y ) = —V{ X) yy<T x <T y + K{ x>y y<rl<Ty 

with Vi x y \ = V > 0, K{ x ,y} = K E R if x, y are nearest neighbors and Vi x y \ = K{ x ,y} = 
otherwise. 

Specifically on the BEG model, there are some rigorous results |11[ [5] concerning the unique- 
ness of the equilibrium state in the regime where the ferromagnetic interaction V is sufficiently 
strong. Our results imply that no phase transition can occur in the BEG model when the crystal 
field is sufficiently large. 

2 Proof of theorem 1 

We will prove the theorem by performing a high temperature polymer expansion on the partition 
function of the system. We will see below that such an expansion, when applied to a spin system 
described by the Hamiltonian (1.1), converges also in the region of low temperatures (large 0) 
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as soon as D > J and the gap between the analytic high temperature phase and the analytic 
low-temperature phase depends on the strength of the crystal field D. When D is sufficiently 
large, this gap disappears. 

We begins by rewriting the partition function (1.5) of the system in the volume A C V as follow. 

ZtSfi) = ^2 e _/3H(CTA) = ^2 e _/3D ^eA CT - e -PJ2{ x ,y}cA v ( x >y>°x<°y) (2.1) 

We then expand the exponential in l.h.s. of (2.1) as follows. 

e -PT,{ x , y }cA V ( x >y> CT x> a 'y) — ^ [ e _ ^S{x,y}cA V(x,y,a x ,a y ) _ ^ + ^ = 

{x,y}cA 

|A| 

= E p{Ri^Ri)--- P(Rs,CTR s ) 

s=l {R 1 ,...,R s }en s (A) 

Here tt s (A) denotes the set of all partitions of A in s non empty subsets, and 

'1 if \R\ = 1 

p(R,a R ) = < 

{ J2 9 eG R U {x , y} eE g [e-^w.,*,) - 1] if \R\ > 2 

with Gr being the set of all connected graphs with vertex set R (E g denotes the edge set of a 
graph g G Gr). Thus (2.1) becomes, 



|A| 

Z A (J3) = Yl e-W*^*Y, E p(Ri,<t Ri )-p(R„<tr.) = 

cr A eS A S=l fli,...,fl s 67T s (A) 

|A| 

= E E m)-p(Rs) 

s=l R u ...,R a en a (A) 



where 



p(R) = p{R)e- pD ^K a * 

Now, due to assumption A. formula (1.2), we have that Y[{ x y y e E g [e~ f3v ^ x ' y ' ax ' ay ^ — 1] = for 
any g £ Gr whenever o x = for some x G R. Hence a straightforward computation gives 



~p{R) 



i + 2Zk=i e ~ (3Dk2 [f \R\ = 1 

£ CTfle £^ DE ^ CT *£ 9eGfl n^^ if \R\ >2 

where Y> R = {o R eT, R : o x / 0, Vie i?}. 
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Define now, for any finite R C V such that \R\ > 2 

\R\ 



C(R) 



l 



-f3Dk 2 



E ■ 



e n i 

g&G R {x,y}<EE g 



-PV(x,y,a x ,cTy) 



1] (2.2) 



Then it is not difficult to check that 



N 



Z A (P) = (1 + 2JV 



-,ODfc 2 \|A 



(2.3) 



k=l 



where Ha(/3) is given by 

Ea(/3) 



1 + ^ n! ^ 



C(Ri) ■ ■ ■ ((Rn) 



(2.4) 



S>1 



|Hj|>2, RiDRj = 

Hence the free energy of the system is given by 

AT 

/AOS) = log(l + 2 ^ e-W) + P k {(3) 
fc=i 



where 



PA(/3) = ^rlogH A ( / 9) 



(2.5) 



The function log(l + 2^^^ e _/3Dfc2 ) is analytic for all > so to check the analyticity of the 
free energy it is sufficient to check the analyticity of the function P\(f3) defined in (2.5). 

Now we just note that S\(j3) is the partition function of a hard-core polymer gas in which 
the polymers R are finite subsets of V with cardinality greater than 1, activity C(R) an d with 
the incompatibility relation being the non-empty intersection. By Fernandez-Procacci criterion 
[6] (see also sec. 3 in [10J ) , the pressure Pa(P) °f this hard-core polymer gas is written in terms 
of an absolutely convergent series bounded uniformly in A if the activities (2.2) satisfy 



mf^-ir^-sup IC(A)I<1 



(2.6) 



n>2 



flCV: xSR 
\R\ = n 



Now, since e ^E^es a 'l < e PO\R\ f or an y aR ^ an( j se tting 



A 



/3 



-f3Dk 2 



(2.7) 



we have 



sup 

x ^ RCV:xeR 
\R\ = n 



C(R)\ < A3 



^sup 



E 

fiCV: zgij 
J?| = n 



E 



e n i 



,-^V(a!,i/,<7a;,cr v ) 



1] 



(2.1 
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We now need to bound the factor 



J^J [ e -pV(x,y,(T x ,<T y ) _ 
geGfl {ie,2/}g£ 9 



(2.9) 



in the r.h.s. of the equation above. Usually factors like (2.9) are bounded using sophisticated 
tools in cluster expansion such as the so called Battle-Brydges-Federbush tree inequality (see e.g 
[2], [16]) or the Brydges-Kennedy-Abdesselam-Rivasseau formula (see e.g [I], [3]). However, for 
the case considered here (bounded spins in a discrete set), it is possible to use a much simpler 
result obtained recently in [10] whose very simple proof is based on the Penrose identity [H] (see 
also [15], [IB], and [6]). Such a result, stated in [10] as Proposition 4.3, in the present notations 
can be rephrased as follows. 

Proposition 2 (Jackson-Procacci-Sokal). Let -P2OO be the set of unordered pairs of distinct 
elements of V and let v : fb(V) — > C : {x,y} 1— > w/j- yi 6e a map u>ii/i i/ie following property: 
there exists a (nonnegative) constant B such that 

sup^l^}! <2B (2.10) 
Then, for any finite R C S with \R\ > 2 the following inequality holds: 
g€G R {x,y}eE g t£T r {x,y}eE T 

where Tr is t/ie sei 0/ all trees with vertex set R and E T denotes the edge set of a tree r G Tr. 

Proof. Choose a root xo in i? so that any tree r G Tr is regarded as rooted in xq- Choose an 
order in i?, i.e. to any 1 in i? associate in a one-to-one way a number in {1, 2, . . . , \R\} (the 
label of x). Then for any tree t eTr with edge set £V define a graph p(r) in Gr with edge set 
-Sp(r) ^ obtained by adding to 2£ r all unordered pairs in connecting each vertex x ^ xq of 
r to all its siblings in r (vertices of the same generation) and to all its "uncles" in r (vertices of 
the previous generation) with label greater than the label of the parent of x. Then the Penrose 
identity is 

H (e-W - 1) = J2 II {e~ V{x ' y} -l)e" E{2 ' z ' }e ^M\^" {z ' z ' } 

g£G R {x,y}€E g reT R {x,y}£E r 

Hence 

I Y, II {e~ v ^,y} - 1) I < Y II (^"^ - J ) e E{z ' z ' }eB ^)\^ k{z ' z ' }l = 

geG R {x,y}£E g t€T r {x,y}&E T 

= V [J (l-e" 1 ^-.!'} 1 ) e E{z ' z ' }eB pM |,;{z,z,}l < X] (l- e -K-.»}l) e S {z ,,' } cfll«{ z , z '}l < 

-£ FI (1 -e" 1 ^,!/}!) e < e ^|R| ^ JJ (!_ e -K^}l) 

t£T r {x,y}£E T reT R {x,y}£E T 

□ 
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Now, for vs XjV \ = (3V(x, y, a x , a y ), by assumption B., formulas (1.4) and (1.3), the pair 
potential V(x,y,a x ,a y ) satisfies (2.10) with B = J. Therefore, proposition 2 leads to the bound 



e n 



-pV(x,y,a x ,(T v ) 



1] 



< e PJ\R\ I] (l-e-^ y) ). (2.11) 

teT R {x,y}eE t 



geG R {x,y}£E g 

Hence, pugging (2.11) in (2.8) and using also that Yla R £E R 1 = (2iV) n , we get 



sup Yl IC(«)I < [2iVV^]«sup Yl E II ( l e 



-PJ(x,yh 



icgV 



RCV: xCR 
\R\ = n 



rcv-.xzr teT R { x ,y}eE t 

it — n 



[2N\ el 3J } n 
(n-1)! x6 v 



ee n (i 



-/^J(xi,x J ■)^ 



(2.12) 



where T n denotes the set of trees with vertex set {1, 2, . . . , n}. It is now easy to check that 



{i,j}e£ t 



with 



(2.13) 



(2.14) 



y^x 



Note that, by (1.4), we have 

h(P, J) = sup V pj(x, y) [ e-M*'** < sup V /3 J(x, y) < 2/3 J 



(2.15) 



y^x 



y^x 



Using finally Cayley formula (i.e. X/teT 1 = n ™ ) f° r t ree counting, we obtain 



71-2 



*ev B ^ B (n-1)! 



|H| = n 



Hence the convergence criterion (2.6) is 



inf (e a - l)" 1 J2 e<in \ 2Nh (P, J)~^e 



a>0 



n=2 



n-1 



n 



n-1 



< 



1 



n\ 2N\ p e< 3J 



(2.16) 



which, by lemma 6.1 in [TO] is satisfied when 

1 



2Nh(J3, J)X^ J 



> mXpe 133 + 3 



i.e. when 



J]P F{I3) > h(J3, J) 



(2.17) 
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with 

1 f-i i o v 



1 (l + 2Ef = i^ Dfc2 ) 2 

2 l8N 2 e-( D - J )P + 37V(1 + 2 J2k=i e"^ 2 )] 



= , n nl „r ^TT (2-18) 



Now it is easy to see there always exist two positive numbers f3\ and (3% (with < /3 2 ) depending 
on D, J and N, such that the inequality (2.17) is in general satisfied either when (3 < (3\, or 
when (3 > (3 2 and D > J. 

Indeed, F(0) = j spj+iPnz ; so that for /3 sufficiently small, say (3 < f3±, the inequality (2.17) 
is always satisfied, because the l.h.s. tends to -F(O) > for (3 — > regardless of the value of .D, 
while in the r.h.s. lmig^o h((3, J) = 0, by (2.15). 

On the other hand, F((3) tends to as (3 — > oo. So, as soon as D > J, the inequality 
(2.17) is surely satisfied for /3 sufficiently large, say /3 > /?2, when the exponential on the l.h.s of 
(2.17) beats the function h((3, J), which, by (2.15), is always below the linear function 2(3 J. 

Finally, let us show that there is a critical value of D c (depending on J) such that, whenever 
D > D c , the inequality (2.17) is always satisfied for all (3 G [0, oo). Indeed, an upper bound for 
D c can be easily found if we are not looking for optimal bounds. As a matter of fact it is easy 
to see that 

F iP) > " o, for all f3 > (2.19) 

Moreover using the bound (2.15) the inequality (2.17) is surely satisfied if 

e {D-j)(3 > ( 12N + 32 N 2 )(3J (2.20) 

which is always satisfied for any (3 > as soon as D > ^1 + 12Af + 32Af2 ^ J 

So we get for the critical value D c of the crystal field the upper bound 

/ 12iV + 32iV 2 \ 
D c < 1 + J (2.21) 



Of course, due to the rough bounds (2.15) and (2.19), this estimate is far from being optimal. In 
particular, the replacement in (2.17) of the function h((3,J) with its upper bound (2.15) seems 
to be a quite crude estimate, especially in the large (3 regime. To find with accuracy the regions 
of validity of the inequality (2.17) so as to improve the upper bound for D Cl one should be able 
compute explicitly function h(f3, J) which depends sensibly on the behavior of J(x, y) and this 
would go beyond the scope of the present note. 

We conclude by noting that improvements on the upper bound (2.21) might also be obtained 
by trying to enlarge the analyticity region in the disordered phase (i.e. when D > J) via a 
genuine low-temperature contour expansion for the system, i.e an expansion around the ground 
state cr = 0. 
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